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ABSTRACT 

In this paper, we present our HARPS radial-velocity data for eight low-activity solar-type stars belonging to the HARPS volume- 
limited sample: HD6718, HD8535, HD 28254, HD 290327, HD43197, HD44219, HD 148156, and HD 156411. Keplerian fits to 
these data reveal the presence of low-mass companions around these targets. With minimum masses ranging from 0.58 to 2.54 Mj up , 
these companions are in the planetary mass domain. The orbital periods of these planets range from slightly less than one to almost 
seven years. The eight orbits presented in this paper exhibit a wide variety of eccentricities: from 0.08 to above 0.8. 

Key words, stars: individual HD6718; HD8535; HD 28254; HD 290327; HD43197; HD 44219; HD 148156; HD 156411 - stars: 
planetary systems - techniques: radial velocities - techniques: spectroscopic 
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^2 , 1. Introduction (IBaranne et al.ll996h . The second fibre can also be used formea- 

. ■ suring the spectrum of the sky background. With more than 70 

CO In this paper, we present new detections from the HARPS planet detections obtaine d so far, including the lightest planet (Gl 58 1 e, 

O search programme. It has been ongoing since 2003. Its de- iMavor et alJl2009h . the HARPS contribution to the excrolanetre- 

. tections are based on radial velocities obtained with the High searc h field is a major one 
■ Accuracy Radial velocity Planet Searcher (HARPS, Mayo r et all 

^ 1 [2003). HARPS is a high-resolution echelle spectrograph mounted Slnce lts beginning, the HARPS planet search programme 

£ i on the ESO-3.6-m telescope in La Silla (Chile). The instrument has been dlvlded lnto several sub-programmes related to spe- 

>< ■ is installed in a vacuum vessel, and it is thermally controlled. clfic sclentlfic questions. One of these sub-programmes is the 

H ! These characteristics result in a very high stability of the intru- HARPS volume-limited programme to which about 20% of the 

mental radial-velocity zero point (measured drifts smaller than HARPS Guaranteed Time Observation (GTO) programme allo- 

lms 1 over one night). HARPS possesses two optical fibres. cated tlme was dedicated. The sample of the HARPS volume- 

The first one collects the target light, whereas the second one llmlted Programme is an extension up to 57.5 pc of the CORALIE 

can be used for recording a reference Thorium- Argon spectrum, ,P lanet search sample (volume-limited sample up to 50 pc, see 

The second spectrum has been used for measuring the residual l Udry et al.|2000|) . It contains about 850 F8-M0 dwarfs. As in the 

drift of the instrument since the time of the latest wavelength case of the CORALIE sample, the maximum distance decreases 

calibration, allowing us to apply a correction to the velocities wlth the B ~ V colour lndex for tar 8 ets wlth late s P ectral ^ 

" (later than K0). Our target list will be published in another paper 

(Lo Curto et al., in prep.). 

* Based on observations made with the HARPS instrument on the 

ESO 3.6-m telescope at the La Silla Observatory (Chile) under the GTO Tar 8 ets ln thls sam P le »« observed at lower signal-to-noise 

programme ID 072.C-0488. ratios (typically 40-50) than the ones in the high-precision pro- 

** Radial-velocity tables are only available in electronic form at the gramme and without using the simultaneous calibration. The 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via typical photon-noise error obtained in this case is about 1-2 m s _1 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?I/A+ A/ depending on the target spectral type and projected rotational ve- 
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locity. This strategy, optimized for the follow-up of a large stellar 
sample, is a good compromise between precision and observing 
efficiency. 

Rather short exposure times (less than 5 minutes) are usually 
enough for reaching the targetted signal-to-noise ratio. An expo- 
sure time of 5 minutes is too short for correctly averaging out in- 
trinsic stellar signals due for example to pulsation or granulation. 
Our velocities can thus be affected by those effects that in some 
cases can reach a peak -to-peak amplitude of several ms" 1 (se e 
e.g. the case of fi Arae, ISantos et alj2004al:lBouchv et al. 2005). 
The amplitudes of these signals depend on the target spectral 
types and on the averaging strategy. Expected radial-velocity er- 
rors due to intrinsic signals are greater for earlier spectral types 
and for evolved targets. 

With such an observing strategy, our programme is mostly 
sensitive to Saturn or Jupiter-like planets but is also sensi- 
tive to lower-mass planets around stars with a low intrinsic 
radial-velocity variability. Also, we generally avoid taking ex- 
posures shorter than 1 minute in order to reduce the impact 
of observing overheads. As a consequence, the brightest tar- 
gets in our sample frequently have photon-noise errors sig- 
nificantly lower than the targetted one. For illustrating these 
facts, we mention th e recent detection s of a 5.5M ffl planet 



around HP 21 5497 ( Lo Curto et alj |2010h , of the 12.6 M ffi planet 
HD 125595 b dSegransan et alJl2010h and o f the 14.4 M e planet 
around BD -08 2823 b dHebrard et al.ll2010h . 

The main aim of this observational effort is a better char- 
acterization of orbital elements distributions of gaseous plan- 
ets. With the increasing size of the sample of known plan- 
ets, a quantitative comparison between the observed sample 
and the predictions of planet for mation models becomes possi- 
ble and reliable (see for example Mordas ini et alj|20 09). With 
the eight planetary companions presented in this paper, the 
total number of planets detected in the volume-limited sam- 
ple now reaches 39 ( s ee iPepe et alj|2004afc [M outou et al. 2005; 

Naef et all l2007t iMoutou et all 120 09b: 
Mordasini et al.l l2010l "TSegransa n et al.l 



Lo Curto et al 



Lo Curto et al. _ 
20 lOt IMoutou et al.ll2010t) 



In this paper, we present our HARPS radial-velocity 
data for 8 solar-type stars hosting low-mass compan- 
ions: HD6718, HD8535, HD28254, HD290327, HD43197, 
HD 44219, HD 148156, and HD 15641 1. This paper is organized 
as follows. We present the main characteristics of the 8 host 
stars in Sect. [2] Section [3] contains the description of the radial- 
velocity sets and orbital solutions obtained for these targets. 
Finally, we provide concluding remarks in Sect. [4] 



2. Characteristics of the host stars 

We list the main characteristics of the nine host stars in Table Q] 
Apparent V-band magnitudes my and colour i ndexes B - V are 
from the HIPPARCOS Catalogue dESAlll997l) . Distances d are 
based on the improved HIPPA RCOS astrometric parallaxes n de- 
rived by Ivan Leeuwenl (|2007). Absolute magnitudes My were 
computed using my and d. 

For each star, we derived the atmospheric parameters T e fi, 
log g, and [Fe/H] by performing a detailed local thermodynam- 
ical equilibrium spectroscopic analysis. A high signal-to-noise 
co-added HARPS spectrum was built for each target. The set of 
Fe I/Fe II equivalent widths used for the anal ysis were measure d 
on the co-added spectra with the ARES code (ISousa et al. 2007). 
The list of spectral lines for th e analysis, as well as more details 
on the method can be found in lSantos et al.l (l2004bl) . 



Stellar luminosities L were computed using HIPPARCOS ap- 
parent magnitudes and parallaxes and the bol ometric correc- 
tion obtained via the r e (j-based calibration in Flower (1996). 
Stellar masses M», and radii R * were derived from r a ff, [Fe/H], 
and My by interpolati ng the Girardi et al.l d2000t) theoretical 
isochrones. Following lFernandes & San tos (2004), a relative un- 
certainty of =s 10% can be assumed for M». We also com- 
puted surface gravities, \ogg, using M„ and R*. These values 
are in fair agreement with the spectroscopic ones: a mean dif- 
ference (spectroscopy-isochrones) of 0.06 dex and an rms of 
the differences equal to 0. 1 1 dex. Spectral types listed in Tabl^U 
are b ased on the s tellar parameters we determined in this pa- 
per. iJenkins et al. (2008) have published iron abundances for 
HD43197 and HD 148156. Their values, [Fe/H] = 0.37 + 0.09 
and 0.22 + 0.15, respectively, are in good agreement with ours. 
The atmo spheric parameters for H P 442 1 9 have also been deter- 
mined by Rob inson et alj d2007l) . They found r eff = 5723 ± 82, 
logg = 4.37 + 0.13 and [Fe/H] = -0.05 ± 0.07. These values 
agree with ours within uncertainties. 

Activity indexes log/? HK were extracted from individual 
HARPS spec t ra usi ng a method similar to the one described in 
ISantos et al.l ([2000). Values indicated in Table[T]are average ac- 
tivity indexes obtained for each target. All the targets exhibit 
a very low level of activity. Our logR HK values for HD43197 
and HD 148156 are in good agreement with those published in 
IJenkins et al.l (12008): -5.12 and -5.08, respectively. We did not 
try to estimate chromospheric ages for our targets because these 
two quantities are no l onger correlated for weakly active stars 
dPace & Pasquinil2"00 4). However, we can infer for most of these 
stars (if not all) a safe age constraint based on their low activity 
indexes and rather high luminosities: age > 3 Gyr. The projected 
rotational velocities v sin i were derived from the widths of the 
HARPS cros s-correlation funct ions (CCF) following the method 
presented in lSantos et al.l d2002l) . 

Finally, HD 28254 is the bright component of a visual bi- 
nary. Accordin g to the Catalogue of Co mponents of Double and 
Multiple stars dDommanget & Nvsl2002l) . the visual companion, 
CCDM J04248-5037B, is located at a separation of 4.3"and at a 
position angle of 254 degrees. With an apparent magnitude of 
my = 13.8, this object is 6.11 mag fainter (i.e. a flux ratio less 
than 0.4%). It was thus clearly outside the detection capabilities 
of HIPPARCOS. Assuming the two stars form a bound system 
and using the distance of HD 28254 (54.7 pc), we computed the 
projected separation between them: 235 AU. We also computed 
the absolute magnitude of the B component: My = 10.11 corre- 
sponding to an M0V-M2V star. Wit h this absolute magn itiude 
and the mass-luminosity relations in iDelfosse et al.l d2000l) . we 
computed a mass of 0.48 M for this object. Its minimum period 
is about 1000 yr. Assuming that this companion is on a circular 
edge-on orbit, the maximal radial-velocity amplitude it would 
induce on the primary star is of the order of 1 . 1 km s _1 . 

3. HARPS radial-velocity data and orbital solutions 

As explained before, stars belonging to the HARPS volume- 
limited sample are generally observed without simultaneous 
thorium-argon reference. Possible instrumental drifts therefore 
remain uncorrected. This has only a limited impact on our re- 
sults because: 1) the instrument is very stable: on average, the 
radial-velocity drifts of HARPS remain below 0.5 m s over one 
night and 2) the targetted photon-noise error for this programme 
is 1-2 m s . We stress again the fact that our observing strategy 
is not optimally designed for averaging out stellar intrinsic sig- 
nals. We list in Table [2] the median exposure times for all the 
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Table 1. Observed and inferred stellar characteristics of the host stars. 



Parameter 


Unit 


HD6718 


HD 8535 


HD 28254 


HD 290327 


5 p. Type 




G5V 


GOV 


G1IV/V 


G8V 


My 




8.45 


7.70 


7.69 


8.99 


B-V 




0.662 


0.553 


0.722 


0.761 


71 


(mas) 


18.23 ±0.76 


19.03 ±0.60 


18.29 ±0.53 


17.65 ±1.57 


d 


(pc) 


54.9!^ 


52 5 +1 ' 7 


54.T\j 


56 7 +5 5 

—4.6 


My 




4.754 


4.097 


4.001 


5.224 


B.C. 




-0.145 


-0.028 


-0.102 


-0.126 


L 


(Lq) 


1.13 


1.86 


2.18 


0.72 


Teg 


(K) 


5746 ± 19 


6136 ±18 


5664 ± 35 


5552 ±21 


logg" 


(cgs) 


4.48 ±0.03 


4.46 ± 0.06 


4.12 ±0.05 


4.42 ± 0.04 




(cgs) 


4.40 


4.35 


4.13 


4.39 


[Fe/H] 




-0.06 ± 0.02 


0.06 ±0.02 


0.36 ±0.03 


-0.11 ±0.02 


M, 


(M ) 


0.96 


1.13 


1.06 


0.90 


R, 


(R Q ) 


1.02 ±0.03 


1.19 ±0.04 


1.48 ±0.06 


1.00 ±0.01 


logR HK 




-4.97 


-4.95 


-5.10 


-4.96 


v sin i 


(kms -1 ) 


1.76 ± 1.0 


1.41 ± 1.0 


2.50 ± 1.0 


1.44 ± 1.0 


Parameter 


Unit 


HD 43197 


HD 44219 


HD 148156 


HD 156411 


S p. Type 




G8V 


G2V 


F8V 


F8IV/V 






8.98 


7.69 


7.69 


6.67 


B-V 




0.817 


0.687 


0.560 


0.614 


71 


(mas) 


17.76 ± 1.22 


19.83 ±0.78 


19.38 ±0.75 


18.25 ±0.49 


d 


(pc) 


56 3 +4 2 

JU - -3.6 


50.4+ 2 ' 


51.6!^ 


54.8+^ 


My 




5.227 


4.177 


4.127 


2.976 


B.C. 




-0.137 


-0.085 


-0.009 


-0.060 


L 


(L ) 


0.73 


1.82 


1.78 


5.38 


T e g 


(K) 


5508 ±46 


5752 ±16 


6308 ± 28 


5900 ± 15 


logg" 


(cgs) 


4.31 ±0.08 


4.21 ±0.03 


4.56 ± 0.07 


4.07 ± 0.04 




(cgs) 


4.42 


4.20 


4.36 


3.87 


[Fe/H] 




0.40 ± 0.04 


0.03 ±0.01 


0.29 ± 0.02 


-0.12±0.01 


M, 


(M ) 


0.96 


1.00 


1.22 


1.25 


R, 


(R ) 


1.00 ±0.01 


1.32 ±0.07 


1.21 ±0.03 


2. 16 ±0.09 


logR HK 




-5.06 


-5.03 


-4.94 


-5.05 


v sin (' 


(kms" 1 ) 


2.18 ± 1.0 


2.22 ± 1.0 


5.70 ± 1.0 


3.30±1.0 



Notes. (r) From our spectroscopic LTE analysis ( 1 Using M, and R, 



targets. This median is below 3 minutes in all cases. We thus ex- 
pect our velocities to be affected at some level by stellar intrinsic 
jitter in particular the velocities of our early type and/or slightly 
evolved targets, namely HD8535 (GOV), HD 28254 (G1IV/V), 
HD 148156 (F8V), and HD 15641 1 (F8IV/V). 

Our radial velocities can be obtained from the extracted 
spectra by cross-correlating them with a numerical template. 
Numerical templates are chosen to minimize the mismatch with 
the target spectral type. Our radial-velocity uncertainties include 
two well-quantified contributions: photon noise and calibration 
error. An additional error of 0.5 ms 1 is added quadratically 
to account for the absence of simultanous calibration and thus 
the possibility of a small uncorrected radial-velocity drift of the 
HARPS zero point. 

All the orbital solutions presented in the following sections 
were obtained with software that includes a genetic algorithm. In 
general, the only a priori information used by this algorithm is 
the type of model we want to use (i.e., 1-Keplerian, 2-Keplerian, 



1-Keplerian plus a linear/quadratic drift, etc.). Once the global 
parameter space minimum is identified, we use it as the starting 
point for a standard Levenberg-Marquardt minimization out of 
which the final solution is obtained. In all the cases, we have ver- 
ified the absence of correlation between the post-fit residuals and 
the measured cross-correlation function bisectors. Finally, un- 
certainties on the fitted parameters are determined using Monte- 
Carlo simulations. 



3.1. HD6718 

The 22 HARPS radial-velocity measurements were gathered for 
HD6718 (HIP5301) between BJD = 2452 991 (December 18, 
2003) and BJD = 2 455 019 (July 6, 2009). These velocities have 
a mean uncertainty of 1.59ms -1 . We have fitted a long-period 
Keplerian orbit to these data. The orbital parameters resulting 
from this fit are listed in Table [2] With these parameters and as- 
suming a mass of 0.96 M for the host star, we computed a mini- 
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Table 2. Orbital solutions for the eight host stars. 



Parameter 


Unit 


HD 6718 


HD8535 


HD 28254 


HD 290327 


P 
T 

e 

y 

OJ 

A, 

Slope 
Curvature 
f(m) 
a\ sin; 


(d) 

(JD" ) 

(kms- 1 ) 

(deg) 

(ms 4 ) 

(ms~' yr" 1 ) 

(ms~' yr~ 2 ) 

(10-"M G ) 

(10- 3 AU) 


2496 ± 176 
357 ±251 
10 +(m 
34.7509 + 0.0013 
286!g 
24.1 ±1.5 

3 56 +071 
5 50 +0 - 45 

JJ -0.24 


1313 + 28 
537 ±92 
15 +009 
2.4588 ±0.0008 
61 +18 

UI -43 

11.8 + 0.8 

n 99+O.O6 
u - iZ, -0.04 

1-41^ 


1116 + 26 

49+43 

U - 8 1 -0.02 

-9.315* ±0.002 
301 +3 
37.3+^ 
2.629 ±0.015 
-0.9290 ±0.00012 

1 99+0.30 
1-i -0.17 
9 oc+0.19 

-0.13 


2443+^ 
1326 +23 ° 

29.559 ±0.003 

26$ + _\t 
41.3 ±2.9 

n 94+0.61 
7 -0.29 


ni2 sin i 

CI 


(M Jup ) 


1 56 +0U 

t SA+0.24 
J.3O_ 015 


0.68!»;S 

7 4^+0.04 
-0.06 


1 16 +010 
1- -0.06 

9 1 1+0.04 

^ -0.05 


9 «;4+o.i7 

o 4^+0.20 
-0.12 


Med(f„ n ) 
(S/N> 
<£rv> 
5 pan 

o~o-c 


(s) 

(ms" 1 ) 
(d) 

(ms" 1 ) 


22 

158 
52 
1.59 
2028 
1.79 


45 
106 
58 
1.57 

2220 
2.49 


32 

60 

51 
1.34 
1989 
2.19 


18 
173 

45 
1.73 
1986 
1.60 


2 




1.50 


2.48 


3.08 


1.22 


Parameter 


Unit 


HD 43197 


HD 44219 


HD 148156 


HD 156411 


P 
T 

e 

y 

(0 

A, 

Slope 
Curvature 

/(«) 
d\ sin i 


(d) 

(JD fl ) 

(kms- 1 ) 

(deg) 

(ms 4 ) 

(ms~ ! yr -1 ) 

(ms~' yr" 2 ) 

(10-"M Q ) 

(10" 3 AU) 


327.8 ± 1.2 
713. 8± 1.7 
0.832g 
72.512+*; 
25 \ +1 
32.4 + '° 3 9 

21 +013 
55+ 010 

U ' J -0.04 


472.3+I 3 , 
585.6^ 2 

Q 61 +0.07 
u ' u -0,09 

-12.0732 c ±0.0026 
147.4+j 2 
19.4 ±3.0 
-0.637 ±0.016 
-0.76390 ±0.000 
17 +005 

'-0.03 

66 +0 06 

u.uu 004 


1027 ± 28 
707 ± 20 

52 +0 04 

"■ J -0.09 
-l./480_ 0009 

35 +17 

JJ -12 

17 5 +L5 

35 +008 

-0.04 

1 41+0.09 


842.2 ± 14.5 
356!« 
0.22 ±0.08 
-38.7512 ±0.0011 

45 +22 

t J -14 

14.0 ±0.8 

9^+0.04 
"""-0.03 

1 06 +00S 


ni2 sin i 


(M Jup ) 


60 +0 12 

"•""-0.04 
Q7+ 01 


n co+o.06 
ios_ 004 


0.85*g 

9 4C+O.04 
^•^ -0.05 


74 +0 05 

/H -0.04 

1 eo+0.03 
i - 58 -0.04 


N 

1 * mes 

Med(f cxp ) 
(S/N> 
(£rv> 
5 pan 

cro-c 


(s) 

(ms- 1 ) 
(d) 

(ms" 1 ) 


33 
171 
48 
1.42 

1943 
1.44 


46 

63 

51 
1.55 
1988 
2.39 


42 

92 
56 
2.24 

2168 

3.69 


50 
90 
73 
1.26 
2231 
2.94 


^red 




1.19 


2.62 


2.89 


5.21 



Notes. Errobars are Monte-Carlo based lcr uncertainties. <6rv) is the radial-velocity uncertainty weighted average. <x -c is the weighted rms of 
the residuals to the fitted orbits. Med(/ exp ) is the median exposure time and (S/N) is the average signal-to-noise ratio at 550 nm. 
w JD = BJD-2 454 000 (W at BJD =2 454 106.6 <f) at BJD =2 454 5 12.6 



mum mass of 1 .56 Mj up for the companion and a semi-major axis 
of 3.56 AU. Residuals to the fitted orbit have a very low disper- 
sion (<r _ c = 1.79 ms -1 ), and a 1 -planet model is thus sufficient 
for explaining the detected signal. We display our data and the 
fitted orbit in panel a of Fig. [TJ 



3.2. HD8535 

A set of 45 HARPS radial velocities was obtained for HD 8535 
(HIP 65 11) between PJD = 2452 850 (July 29, 2003) and 
BJD = 2455 070 (August 26, 2009). The mean radial-velocity 
uncertainty of these measurements is 1.57ms~'. The presence 
of a long-period signal is obvious (see panel b of Fig . [TJ . The re- 
sults of the Keplerian fit to our data are presented in Table|2] The 
dispersion of the residuals to the fitted orbit are somewhat larger 
than expected (y? , = 2.48). We have tried to fit orbits including 
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BDJ - 2 450 000 [d] 

HD 148156 


5000 
HARPS 
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r I 1 
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3000 3500 4000 4500 
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HD 156411 
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Fig. 1. HARPS radial-velocity data for the eight stars presented in this paper. Temporal velocities (filled dots) and fitted orbital 
solutions (solid lines) are displayed in the upper panels. Residuals to the fitted orbits are displayed in the bottom panels. For 
HD 28254 and HD44219, the dashed lines in the upper panels show additional quadratic radial-velocity drifts that were included to 
the orbital fits. For all the other targets, the upper panels dashed lines represent the systemic velocity y. 



additional signals (linear or quadratic trend; 2-Keplerian model) 
but we were unable to significantly improve the quality of the 
solution. We have also verified the absence of any correlation 
between the residuals and the CCF bisectors. 



Assuming a primary mass of 1 . 1 3 M Q , we computed the min- 
imum mass of the companion and the semi-major axis of its or- 
bit: m.2sini = 0.68 Mj up and a = 2.45 AU. The fitted orbit and the 
fit residuals are displayed in panel b of Fig.Q] The phase-folded 
orbit and velocities are displayed in panel a of Fig. [2] 



3.3. HD 28254 

We collected 32 HARPS radial velocities for HD 28254 
(HIP 20606) between BJD = 2452941 (October 28, 2003) and 
BJD = 2454931 (April 9, 2009). The average velocity uncer- 
tainty of these measurements is 1.34 ms -1 . As indicated in 
Sect. [2] HD 28254 is the main component of a visual binary 
system. The presence of a nearby visual companion can induce 
radial-velocity errors via the seeing-depend ent contaminat ion of 
the spectrum of the main component. iPepe et alj J2004b) have 
simulated this effect for HARPS and using the La Silla seeing 
statisitics. We have performed the same simulations using the 
characteristics of the HD 28254 system (p-4.3", Amy =6.11). 
The expected effect is negligible in this case («0.1ms-'). 
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A long-period signal is visible on the data (see panel c of 
Fig. Q]). From the visual inspection of these data, it is also ob- 
vious that the Keplerian signal is highly eccentric. We first at- 
tempted to fit a single orbit to our velocities. With <?>0.99, 
the obtained orbit was not physical so we rejected it. After 
this initial unsuccessful attempt, we tried to fit orbits with ad- 
ditional signals. The following two models did not give satis- 
factory results: the 1-Keplerian plus linear drift model and the 
2-Keplerian one. The best fit we obtained and finally adopted 
is based on a 1-Keplerian plus quadratic drift model. Despite the 
poor coverage of the second radial-velocity maximum, hence the 
poorly constrained velocity semi-amplitude, the planetary nature 
of HD 23254 b is well established. The orbital period is very well 
constrained. Solutions with significantly larger amplitudes are 
possible, but they also have much larger eccentricities. The ef- 
fect of the increasing amplitude on the minimum mass of the 
companion is largely compensated for by the increasing eccen- 
tricity as the mass function scales with K^(l - e 2 ) 3 ^ 2 . 

More data are necessary for better constraining the shape of 
the orbit of HD 28254 b. In particular, a better coverage of the 
periastron passage will allow us to constrain the radial-velocity 
semi-amplitude and orbital eccentricity. These future observa- 
tions will also help us characterizing the secondary signal bet- 
ter. The orbital solution we present here is thus a preliminary 
one. The observed additional quadratic signal remains compati- 
ble with the signal that the visual companion described in Sect. [2] 
could induce. 

Using our fitted parameters, we computed the minimum 
mass of the planetary companion and its orbital semi-major axis: 
ni2sini = 1 . 16 Mj up and a = 2. 15 AU. Our radial-velocity data and 
the fitted orbit are displayed in panel c of Fig. Q] The phase- 
folded orbit and velocities are displayed in panel b of Fig. [2] 

3.4. HD 290327 

Between BJD = 2452945 (November 1, 2003) and 
BJD = 2454931 (April 9, 2009), we have observed HD 290327 



(HIP 25 191) 18 times with HARPS. The obtained radial veloci- 
ties have a mean uncertainty of 1 .73 m s -1 . A long -period signal 
is well visible on the data. We list in Table [2] the parameters 
resulting from a Keplerian fit to the data. With these parameters 
and using the primary mass of Table Q] (0.90 M ), we computed 
the minimum mass of the companion that induces the detected 
radial-velocity signal and the semi-major axis of its orbit: 
m.2sim : = 2.54 Mj up and a = 3.43AU. The dispersion of the 
residuals to the fitted orbit is very low: 1.6 ms -1 . Our velocity 
data and the fitted orbit are displayed in panel d of Fig.Q] 

3.5. HD43197 

We have gathered a set of 33 radial velocities for HD43197 
(HIP 29550). These data were obtained between BJD = 2452 989 
(December 15, 2003) and BID = 2 454 932 (April 10, 2009), and 
they have a mean velocity uncertainty of 1.42ms -1 . A visual 
inspection of the data displayed in panel e of Fig. [T]shows a pe- 
riodic highly eccentric signal. The parameters resulting from a 
Keplerian fit to the data are listed in Table [2] The fitted eccen- 
tricity is indeed very high and very well constrained in spite of 
our poor coverage of the radial-velocity minimum. The rms of 
the residuals to the fitted orbit is very low: 1.42 ms -1 . This 1- 
planet fit is thus satisfactory. With a primary mass of 0.96 M Q , 
the computed minimum mass of the companion is 0.6 Mj up . The 
orbital semi-major axis is 0.92 AU. Our data, the fitted orbit, and 
the residuals are displayed in panel e of Fig. Q] Panel c of Fig. [2] 
shows the phase-folded orbit. 

3.6. HD44219 

The HARPS data set we have in hand for HD44219 
(HIP 301 14) consists of 46 radial velocities measured between 
BID = 2 452 944 (October 3 1 , 2003) and BID = 2 454 932 (April 
10, 2009). These velocities have an average uncertainty of 
1.55ms -1 . We first fitted a 1 -Keplerian model to these data. The 
residuals to this fitted orbit were a bit large (<x _ c = 2.95 m s -1 ). 
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We then tried to include an additional signal in the fitted 
model. Adding a linear drift improved the quality of the fit in 
a marginally significant way (xLa decreasing from 3.79 to 3.30; 
f-test probability = 71%). The improvement resulting from in- 
cluding a quadratic trend in the fitted model was more signifi- 
cant: xLa = 2.62 and f-test probability equals 90%. We also tried 
to fit a 2-planet model but we were unable to find any convincing 
solution. We finally decided to adopt the Keplerian + quadratic 
drift model. The parameters resulting from this fit are listed in 
Table[2] The orbital parameters obtained with our initial 1 -planet 
fit did not differ much from the ones in the adopted solution 
(P =447 d; e = 0.57; K x = 17 ms- 1 and m 2 sin i = 0.52 M Jup ). 

The dispersion of the residuals to the finally selected solution 
is still a bit large (<xo-c = 2.39 m s -1 ), but this is probably due to 
the still not very well constrained additional signal. With this 
solution and assuming M, = 1 M , the computed minimum mass 
for HD44219b is 0.58Mj up and its orbital semi-major axis is 
1.19 AU. We display the fitted model and the temporal residuals 
in panel f of Fig. Q] The phase-folded orbit and velocities are 
displayed in panel d of Fig. [2] 

3.7. HD 1481 56 

We have collected 42 radial velocities with HARPS for 
HD 148156 (HIP80680) between BJD = 2452852 (July 31, 
2003) and BJD = 2455 020 (July 7, 2009). The mean radial- 
velocity uncertainty of this data set is 2.24 ms -1 . A long-period 
signal is visible in the data displayed in panel g of Fig. Q] The 
results of the Keplerian fit performed on our velocity set are 
presented in Table [2] The rms of the residuals is a bit large 
(o"o-c = 3.69 ms" 1 ;^ = 2.89). We were unable to significantly 
improve the fit by adding extra signals. No obvious correlation 
between CCF bisectors and residuals could be detected. The ori- 
gin of these abnormal residuals remains unknown, but they may 
be due to intrinsic stellar jitter. 

The minimum mass of HD 148156b is 0.85 Mj up and the 
semi-major axis of its orbit is 2.45 AU (assuming M* = 1 .22 M ). 
The fitted orbit and the fit residuals are displayed in panel g of 
Fig. Q] The phase-folded orbit and velocities are displayed in 
panel e of Fig. [2] 

3.8. HD 156411 

We obtained 50 HARPS radial-velocity measurements of 
HD 156411 (HIP 84787). These measurements were taken be- 
tween BJD = 2452 850 (July 29, 2003) and BJD = 2455081 
(September 6, 2009) and they have a mean velocity uncer- 
tainty of 1.26 ms -1 . A long-period signal is visible on the data 
displayed in panel h of Fig. Q] We fitted a Keplerian orbit 
to these data. The resulting parameters are listed in Table [2] 
The dispersion of the residuals to the fitted orbit is too large 
(o"o-c = 2.94ms -1 ), and the fit quality quite poor (xL* = 5.21). 

We tried to fit models including additional signals. First, we 
tried to add linear and quadratic drifts. In both cases, the fit qual- 
ity improved a bit (xLa =4.80 and 4.70, respectively) but the cor- 
responding f-test probabilities indicated that the achieved im- 
provements were only marginally significant (at the 1 cr level). 
A 2-planet fit was also tried, but we were unable to find any 
convincing solution. We also verified the absence of correla- 
tion between CCF bisectors and residuals. HD 15641 1 is slightly 
evolved so we speculate that intrinsic stellar signals are respon- 
sible for the abnormal residuals. 



For this paper, we finally decided to adopt our initial 1- 
planet orbital solution. With this solution and assuming a mass 
of 1.25 M Q for the host star, we computed a minimum mass 
of 0.74 Mj up and an orbital semi-major axis of 1.88 AU for 
HD 15641 1 b. The fitted orbit and the temporal residuals are dis- 
played in panel h of Fig. [1] The phase-folded orbit and velocities 
are displayed in panel f of Fig. [2] 

4. Conclusions 

In this paper, we have presented our HARPS radial-velocity data 
for 8 low-activity solar-type stars orbited by at least one long- 
period planetary companion. HD6718 has a 1.56 Mj up compan- 
ion on a 2496-d low-eccentricity orbit (e = 0.1). The compan- 
ion orbiting HD8535 is a 0.68 Mj up planet with a period of 
1 3 1 3 day s and an eccentricity of 0. 1 5 . HD 28254 b is a 1 . 1 6 M Jup 
planet on a 1116-d highly eccentric orbit (e = 0.81). The com- 
panion to HD 290327 has a minimum mass of 2.54 Mj up . Its 
2442-d orbit is nearly circular (e = 0.08). As for HD 28254b, 
the 0.6 Mj up companion to HD43197 is on a highly eccentric or- 
bit: e = 0.83. Its orbital period is 327 d. HD44219 is orbited by a 
0.58 M Jup planet on a 472-d 0.61 eccentricity orbit. HD 148156b 
is a 0.85 Mj up planet orbiting its parent star on a 1027-d eccen- 
tric orbit (<? = 0.52) and the companion to HD 156411 is a 0.74 
Jupiter-mass planet with a mildly eccentric (e = 0.22) 842-d or- 
bit. 

We find some evidence of additional longer-period compan- 
ions around HD 28254 and HD 44219. In both cases, we are un- 
able to fit full double-Keplerian orbits because of our short ob- 
serving time span. Future observations of these two systems will 
probably allow us to better characterize the companions respon- 
sible for these signals. 

Finally, with the addition of HD 28254 b and HD 43 197 b, 
the population of exoplanets with extremely high eccentric- 
ities now amounts t o 6. The other pla nets with e > 0.8 are 
HD4113 (e = 0.90 3, iTamuz et all 12008b . 1 56846 (e = 0.8472 , 
iTamuz et al.l l2008h . HD 207 82 b (e = 0.92 iJones et al.l 12006). 
and HD 80606 b (e = 0.934, iMoutou et al.ll2009al) . Fiveof these 
stars hosting hi ghly eccent r ic plan ets are members of wide visual 
binaries (see ITamuz et all (120081) for HD411 3 and HD 156846 
lDesidera & Barbierild2007h for HD 20792, and lNaef et all (|2001|) 
for HD 80606).Only HD43197 is not known to be in a multiple 
stellar system. This large number of eccentric-planet hosts in bi- 
naries tend to favour scenarios in which the extreme orbital ec- 
centricity resul ts from th e influence of a third b ody via Kozai 
oscillations dKozail 11962b or Kozai migration dWu & Murravl 
I20031 
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